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Abstract 

We present new regular static isolated cylindrically symmetric solutions for 
the Ginzburg-Landau model which have finite Gibbs free energy. These con- 
figurations (which we call the flux tube and type B solutions) are energetically 
favorable in the interval of the external magnetic fields between the thermo- 
dynamic critical value H c and the upper critical field H C2 which indicates that 
they are important new elements of the mixed state describing a transition 

from vortices to the normal state. 
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The existence of vortices was predicted by Abrikosov JT| within the framework of the 
phenomenological Ginzburg-Landau (GL) theory of superconductivity. Nielsen and Olesen 
0, noticing the mathematical equivalence of the GL model to the Abelian theory of coupled 
gauge and Higgs fields, demonstrated the relevance of the vortex (string-like) structures in 
the high energy particle physics. An isolated vortex is a static regular cylindrically symmetric 
solution of the classical second-order GL equations with the finite energy (line density), see 
earlier numerical results in ||. Our purpose is to present new exact regular solutions with 
the finite Gibbs free energy. We consider isolated structures and thus confine ourselves to 
the case of cylindrical symmetry. 

The Lagrangian of the model is 



where $ is a complex scalar field, F^ v = d^A u — d u A^ is the electromagnetic field strength, 
Dft = <9 M + iAp, while A and /x are the coupling constants. Since /x has the dimension of 
an inverse length, it is convenient to introduce for the cylindrical system (p, 6, z) a new 



use the cylindrically symmetric ansatz A 2 = f(r), $ = (/x/ vA) ( f( r ), where /, cp are two real 
functions. Notice that we do not include the phase factor e m0 for the scalar field (cf. ||). 
The field ip is defined up to a gauge transformation and we find it more convenient to work 
in such a real gauge. The magnetic field has only one component in ^-direction, given by 
H = (/x 2 /A) h, where we denote the dimensionless magnetic field by h — f'/r. 
With the above assumptions, the equations of motion reduce to 




(1) 



(dimensionless) radial coordinate r via p = (\/A//x)r. Looking for static configurations, we 



r 2 f" - rf> = rVf, 



(2) 




(3) 



From ([]]), the energy per unit length for static cylindrical configurations reads 




o 
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Fields which satisfy ©-© also produce an extremum of £. The latter has an absolute 
minimum £ = for the trivial solution {jp = 1, / = 0) which describes the Meissner state 
with superconducting order in all points of a sample and the magnetic field absent inside. 
The vortex |1|-|J presents an example of a nontrivial finite £ solution with the magnetic field 
penetrating along a thin non-superconducting core around r = 0. 

Although £ > for the vortex, one can prove that it becomes energetically more prefer- 
able, above the first critical field H Cl , than the trivial Meissner state by inspecting the Gibbs 
free energy (per unit length) which is defined by Q = £ — / d 2 x (HH cxt ). We will assume 
that the external constant magnetic field is also directed along the z-axis, and use the di- 
mensionless variable defined, as above for H, by H cxt = (fi 2 / X)h ext . For the understanding 
of a transition to the normal state, of particular interest is the functional 

,2 



oc 



AQ = I (In- 



to 



(h - h cxt ) 2 



+ M 2 + ^ + ^ 4 -2^)j, (5) 

which describes the difference Q — Q n h between the Gibbs free energy of a particular config- 
uration and that of the normal state (<p = 0, h = h cxt ). It is important to notice that (jS)) 
has the same equations of extremals as the energy functional (|4]), namely (0)-(^[). However, 
unlike the strictly positive £, the functional (|5]) can have any sign. 

Let us study the regular solutions of the GL equations (@)-(@). At the symmetry axis, 
we find two types of regularity conditions for / and <p: 

(A) Potential / is non-zero while scalar field if vanishes at the origin, 

f = N + ar 2 + -b 2 r 2N+2 + 0(r 2N+4 ), (6) 

J AN(N + 1) V ; ' y J 

(p = br N (l + , - (a - —) r 2 + 0(r 4 )) , (7) 
^ V 2(iV+ 1) V 2NJ V J J ' K J 

where iV = ±1, ±2, ... and parameters a, b are arbitrary. 

(B) Potential / vanishes while scalar field ip is nontrivial at the origin, 

/ = ar 2 (l + b 2 r 2 /8 + 0(r 4 )) , (8) 
3 



b (l + A(6 2 -l)r 2 /4 + 0(r 4 )), (9) 



with some parameters a, b. 

The vortex solution satisfies (A) and is distinguished by the well known asymptotic 
conditions at infinity (</? — > 1, / — > 0). Analysing the (B) case, one can prove that no finite 
energy solutions for type B regularity conditions exist. However, it is possible to find finite 
Gibbs free energy regular type B solutions which satisfy at infinity 



K r ) = -f\ r ) 

r 



— ► h cxt , (10) 

■ oo 

^(r)| r ._oo — ^ 0. (11) 



The simplest solutions with (p having no nodes are shown in Fig. [I] (in all figures we use 
k := vA)- There exist also more nontrivial solutions with nodes. A principal difference of 
these solutions (which we will call "type B" solutions for brevity) from the vortex solutions 
with a quantized flux lies in the fact that the magnetic field is now asymptotically constant, 
and hence the total magnetic flux is infinite. However a reasonable replacement of the flux 
is provided by 

oo 

M := J pdpd6{H - H cxt ) =2n J drr(h - h cxt ). (12) 

o 

Defined formally as a difference of fluxes, this variable is usually interpreted as a magnetiza- 
tion per unit volume |4j] . Unlike the quantized flux for the vortices, M can have an arbitrary 
value for the type B solutions. 

Numerical integration reveals the following interesting properties of the nodeless type B 
solutions: For A > 1/2 they all have AQ < 0, for A = 1/2 always AQ = 0, and AQ > for 
A < 1/2. Hence, for an ideal type II superconductor the type B configuration is energetically 
more preferable than the normal state. However, the type B solutions do not exist for every 
value of the external field h ext . For a given A > 1/2, one finds h cxt < A for all solutions, and 
thus the upper limit is exactly equal to the second critical field h C2 = A. 

However, this is not the end of the story. Besides the type B configurations, there are 
other nontrivial solutions with finite values of the Gibbs functional (f|). Namely, let us 



take the type A regularity conditions (@)-(0) at the origin, while at infinity we consider the 
asymptotics (|iO|)-(|TT|), which suggests a possible physical interpretation of such solutions 
as emerging from a "gluing" of a vortex with a type B configuration. These solutions are 
shown in Fig. 0. Their form explains why we call them flux tubes 0: There is a non- 
superconducting core filled by the magnetic field (looking like a vortex), surrounded by a 
superconducting tube (almost completely free of the magnetic field). Outside such a tube 
the sample quickly reduces to a normal state (if = 0) with a constant external field in it. 

Each family of the flux tube (FT) solutions has two branches. For the nodeless A = 1 
case, we find that one of these branches is characterized by the positive Gibbs free energy, 
and another has negative Gibbs free energy. However, for A > 9/2 both branches describe 
configurations with AQ < 0. We find it convenient to depict these branches in the form of 
magnetization curves, see Fig. |3|. For A < 1/2, both branches have AQ > 0. 

Similarly to the type B configurations, we find the flux tube solutions with nodes. For 
A = 1, all the FT solutions with one node have AQ > 0. However, with increasing A this 



changes. The crucial point is the position of the thermodynamic critical value h c = yA/2 
relative to the "limiting values" of h ext at which a FT family has vanishing magnetization 
and Gibbs free energy. Numerically, we discovered that these limiting points are located, 
on the h ext axis, at A, |A, |A, |A, . . .. Below we explain these values with the help of the 
linearization analysis. 

Since both, the type B and the FT solutions, have a well defined finite Gibbs free energy, 
it is natural to compare them. We find that always A(?ft < AQ-q. For illustration, we 
display AQ as a function of h cxt for A = (2.25) 2 in Fig. |4] which shows that the flux tube 
configurations are energetically more preferable. This holds true for any A. 

The best way to understand the structure of the type B and FT solutions in the limit of 
vanishing magnetization (M — > 0) and Gibbs free energy (AC? — > 0) is to study the linearized 
GL equations. The linearization of the system (@)-(0) leads to the Schrodinger-type equation 
for if with the potential of a circular oscillator. Regular solutions exist only when 
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where n = 0,1,2,... and s = ±1. Corresponding eigenfunctions ip Uj N are given in terms 
of the Laguerre polynomials, with n equal to the number of their zeros (nodes). Let us 
denote hk ■= X/ {2k + 1), k = 0, 1, 2, ... . For n = and sN = —\N\, one finds that the 
maximal eigenvalue of (|13j) is h = sX = sh , which is precisely the second critical field 
h C2 = X. At the first sight, this fact may seem to be well known in the literature on the 
type II superconductors, cf. [|l],f|]. However, the peculiar point is the conclusion that the 
vortices apparently have nothing to do with the determination of h C2 . As we see, the latter 
is determined by the existence condition for the nodeless type B and FT solutions which 
are energetically more preferable than the normal state. It is clear that the eigenfunctions 
y^o.jv = r ' iV ' exp(— \r 2 ) describe the linearized regular type B (for iV = 0) and FT (for N ^ 0) 
solutions. It seems worthwhile to notice also that the starting point for the well known (see, 
e.g., construction of a vortex lattice near h C2 is not the vortex, as one might think, but 

a linearized type B solution. 

The physical meaning of the lower eigenvalues (|13j) is now also clear: these define the 
limiting values for the external magnetic field at which the exact type B and FT solutions 
become "linearizable" and thus disappear with vanishing magnetization and Gibbs free 
energy. This confirms our numerical observations: For instance, both the nodeless FT and 
the type B configurations have the limiting magnetic field values hi = A/3 and ho = h C2 , 
while the FT and the type B with one node "live" between /i 2 = A/5 and hi — A/3. In 
general, the limiting points for the solutions with k nodes are h^+i and hk. As we already 
mentioned, for a given A the position of the thermodynamic critical value h c = JX/2 relative 
to the interval [h^+i, h^} is important. If h c belongs to this interval, then hk+i < h cxt < h^ for 
all solutions in this family. However when h c does not belong to this interval, then [hk+i, h^} 
is extended up to h c . In particular, we find that for every FT family, the magnetization 
diverges when /i ext — > h c . Moreover, one can show that the Gibbs free energy for a vortex 
and for a flux tube solution are equal for /i ext = h c , 0. 
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One may further weaken the conditions at infinity, requiring / — > for the potential, but 
for the scalar field simply demanding regularity and finiteness. This gives another class of 
regular configurations which we call the oscillating solutions (OS). For large r we can neglect 
the nonlinear terms in ([$]), and get for the scalar field the linearized equation (p" + (p'/r + 
\(p = 0. Its solution is the Bessel function ip = J (\/~\r) which is confirmed by the direct 
numerical integration. Such OS seem to be relevant to observations made in ||[7|]. In our 
opinion these new solutions are unstable configurations preceding to the completely formed 
Abrikosov vortices. Notice, that an oscillating character of such solutions may resemble an 
"intermediate" superconductor state with coexisting normal and superconducting regions. 
But a big difference is that the magnetic field penetrates only at the center, and the magnetic 
flux is quantized precisely as for the vortices. Oscillating solutions exist both for the type A 
and B conditions, in the latter case the magnetic field is completely zero. A curious feature 
is that for a finite sample in absence of an external magnetic field an oscillating state is 
energetically more preferable than a purely normal state. 

In each class of solutions the decisive role is played by the values of the parameters (a, b) 
appearing in the regularity conditions (||)-(0) and Numerically, all the reported new 

solutions were obtained after a "fine tuning" of these parameters. Let us draw a kind of a 
"phase diagram" on the (a, b) plane showing explicitly the domains of existence for different 
solutions. Since the vortices, flux tubes, and oscillating solutions all satisfy the type A 
conditions ©-(0), we can display them on the same (a, 6)-plane, see Fig. |] for the case of 
A = 1. The big dot shows the "position" of a vortex, while each curve represents a complete 
family of solutions (a point on a curve gives (a, b) for a particular solution). The curves hit 
the ct-axis (b = 0) at one of the possible limiting values: a = ±|/ifc, k — 0, 1, — Curves for 
FT with increasing number of nodes are "concentrating" around the oscillating curve which 
seems to indicate that the oscillating solutions are unstable and a small perturbation may 
cause their decay into a nearby flux tube with a finite number of nodes. Moving along any 
FT curve away from the a-axis, one hits the vortex dot, where magnetization diverges. 

In this paper we have presented new numerical solutions of the cylindrically symmetric 
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GL equations. Noticing that the latter are common extremals for the energy (|4]) and the 
Gibbs free energy (|5]) functionals, we describe the regular solutions with finite Gibbs free 
energy. We show that h C2 is not the magnetic field below which the vortex becomes ener- 
getically more preferable than the normal state, but it is the field at which the formation 
of the type B and flux tube structures starts, since A£ FTB < for h ext < h C2 . Our results 
show that the flux tube solutions without node remain energetically most preferable from 
h C2 down to h c , after which the vortices become energetically dominant. It is worthwhile to 
stress, that our results do not contradict the previous knowledge about the mixed state in 
the type II superconductors. Notice that, after all, one can interpret a flux tube solution 
as a vortex "surrounded" by a type B configuration. In a certain sense, one can speak of 
a "fine structure" for the mixed state: Different configurations (vortices, type B and FT 
solutions with k nodes, k = 0, 1, . . .) can exist for any external field h ext between h ci and h C2 
(cf. Table |). 

We are grateful to Friedrich W. Hehl and Eckehard W. Mielke for extremely helpful advice 
and discussions. The work of Y.N. O. was supported by the Deutsche Forschungsgemeinschaft 
(Bonn) grant He 528/17-1, and for F.E.S. by the European Union. 
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FIGURES 

FIG. 1. Type B solutions without node for k = 1.0: The scalar field (p and the magnetic field 
h with initial values <p(0) = 0.9, 0.6, 0.3 and /(0) = 0.0. 

FIG. 2. Flux tube solutions without node for k = 0.5, 1.0, 1.5 with the fixed value of magneti- 
zation. 

FIG. 3. Magnetization curves for nodeless FT with k = 0.5, 1.0,1.5. For each k, the corre- 
sponding limiting values hi, hi, h c , and h C2 are drawn. The magnetization diverges at h c . 

FIG. 4. AG against the external magnetic field for the flux tubes and the type B solutions with 
k = 2.25. The flux tubes have lower AG and thus are energetically more preferable. 

FIG. 5. (a, 6)-diagram for k = 1.0. The big dot describes the corresponding vortex solution. 
The drawn, broken and dotted lines represent the nodeless flux tubes, flux tubes with one node, 
and oscillating solutions, respectively. 
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FIGURE 1: 



. 2 




1 

- 


1 1 1 


i | i | 


- ,v 




h " 






K=1.0 




\ \ cp 




1 




-1 , 1 



FIGURE 2: 





i | i | 


i i i i i i 


I i I 


1 .0 




M/(2tt) = 






-6.0 : 


.5 














h 




if / '' 

y / / * 











K = 


0.5 : 


-.5 




K = 


1.0 : 






K = 


1.5 : 


1 n 








iii" 



2 4 6 8 10 12 14 

r 



FIGURE 3: 




-.5 .5 1.0 1.5 2.0 2.5 3.0 
h e xt 



11 




12 



TABLES 



TABLE I. Existence domains of different GL solutions: fine structure of the mixed state? 

region energetically preferable 

h C3 < /lext normal 

h C2 < h cxt < h C3 surface superconductivity 

hi < h cxt < h C2 0-node FT 

h 2 < h cxt < hi 0-node FT, 1-node FT 

H < h cxt < h 2 0-node FT, 1-node FT, 2-node FT 

h Cl < h cxt < h c pure vortices 

< h cxt < h ci Meissner 
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